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ABSTRACT 

The  calculation  of  vibrations,  and  in  particular,  resonaiKes  from  bounded  elastic  shells  can  be  quite  tedious  and  time 
consuming  when  using  the  exact  elastodynamic  equations.  Thus,  a  popular  approach  has  been  us  employ  vanous 
dynamic  assumptions  about  the  motion  of  the  shell  surface  when  subjected  to  disturbances.  This  can  be  done  using 
variational  considerations  in  which  energy  is  minimized  when  vanous  constraints  art  imposed  We  exploit  the 
technique  using  various  assumptions  which  give  rise  to  several  shell  theories.  We  can  use  the  resulung  expressions 
to  calculate  resonances  over  a  frequency  range  and  compare  them  with  the  exact  results.  We  may  then  rank  the 
various  approximations  in  order  of  (heir  agreement  to  the  exact  results.  Limitations  of  each  of  the  methods  can  then 
be  outlined  as  well  as  those  of  shell  methods  in  general. 

KEYWORDS 

Shell  theory:  elastic:  variational  methods:  nontorsional 
INTRODUCTION 

The  standard  assumptions  used  in  shell  theory  were  formulated  by  A.  E.  H.  Love  (Love.  1944)  and  are  as  follows: 
(I)  The  thickness  of  a  shell  is  small  compared  with  the  smallest  radius  of  curvature  of  the  shell;  (2)  The 
displacement  is  small  in  comparison  with  the  shell  thickness:  (3)  The  transverse  normal  stress  acung  on  planes 
parallel  to  the  shell  middle  surface  is  negligible;  (4)  Fibers  of  the  shell  normal  ua  the  middle  surface  remain  so  after 
deformation  and  are  themselves  not  subject  to  elongation.  We  use  these  assumptions  in  the  development  of  a  shell 
theory  for  an  elastic  sphencaJ  shell  in  the  spirit  of  Timoshenko-Mindlin  plate  theory 

DERIVATION  OF  EQUATIONS  OF  MOTION 

In  spherical  shelb  membrane  stresses  (proportional  to  )  predominate  over  flexural  stresses  (proportional  to  /?’ ) 
where 


I  We  differ  from  the  standard  derivation  for  the  sphere  (lunger  and  Feit.  1986)  by  retaining  all  terms  of  order  in 

both  the  kinetic  and  potential  energy  parts  of  the  Lagrangian.  We  note  that  this  level  of  approximation  will  allow 
I  us  to  include  the  effects  of  rotary  inertia  and  shear  distortion  in  our  shell  theory.  We  begin  our  derivation  by 

I  considering  a  u,v,w  axis  system  on  the  middle  surface  of  a  spherical  shell  of  radius  a  (measured  to  mid-shcU)  with 

thickness  h ,  as  shown  in  Fig.  1 . 
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Fig.  i.  Spherical  shell  showing  ihe  coortiinaie  system  user). 


Lapanpan  Variational  Analysis 

Thus  the  new  Lagrangian  (which  is  equivalent  lo  a  Tunoshenko-MinrUin  theory  as  applied  to  a  spherical  shell)  is 

L  =  T-V*W.  (2) 


where  the  kinetic  energy  is 


(3) 


with  the  surface  displacements  taken  to  be  linear  as  in  Timoshenko- Mindlin  plate  theory; 

u,  =(l  +  -)4---^.  (4) 

a  a  do 

w,=w.  (5) 

There  is  no  movement  in  the  v-dirtcuon  since  the  sound  field  can  be  assumed,  without  loss  of  generality,  as 
torsionless.  By  subsuiuuon.  the  kinetic  energy  is 

or.  simplifying. 

r=  rrp.ha'IJ  ((1.8/3*  +6^‘  +  l)u*  -(3.6^*  +  +  (1.8^*  ■'-(^'  +  l)wM sin  fide.  (7) 

which  to  order  is 
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T  -  *  (I  ♦  ^')'‘Msin  Bd6- 

^0  ^0  Afk 


de 


(8) 


-'■  i 


In  a  similar  fashion  the  potential  ertergy  is 

-*■  edxdm. 


which  by  substitutior  becomes 


f*/l 

r  £  If 

J-4/1 

.1-  v’  (X-Kj)'i, 

x.du  X  , 

((]  +_) - -  +  H.' 

a  dd  a  30^ 


,  />r/.  X  0W  x^du  xdw 

+  (cot  [(1  +  -)tt - — ]  +  h/H(I  +  -)— - ^  +  w] 

a  a  00  a  00  a  00 


(x  +  aY  sin  0dxd0p, 


or  finally. 

y  =  j  f *  ((w  +  +  (w  +  ucot  ef  +  2  vfw  +  +  ucotfl)  + 

I  -  \r  00  00  00  06 

.2vcot8(u-:gK^-|^)l)sm<W9. 


00 


where  the  nonvanishing  components  of  the  strain  are 

e, 


and 


with  nonzero  stress  components  are 


ox) 


where  £  is  Young's  modulus. 


If  0u  'l  X  f  Ai  ^'w'l 

1  ,  ,  X  f  <?w  'I 

=  — (coteu+  w)+— rcotflu-  —  L 
a  a  \  06  J 


I- V 


■(£»  +  ve^). 


(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


Finally,  the  wort  done  by  the  surrounding  fluid  on  the  sphere  is 

W  ^  Imo’  f  Vwsin  006, 

*0 

where  p,  is  the  pressure  at  the  surface  of  the  shell. 

Lagrangian  Densiw  and  Eauarions  of  Motion 

Since  the  integration  along  the  polar  angle  is  intrinsic  to  the  problem,  the  solution  must  be  found  using  a 
Lagrangian  density: 

L,  =  np,ha^[(,l  +  6P^)u^  -  6P\~  -t  (1  +  ^*)w^]sin9  -  ♦(w  +  ucoifi)’ 


00  00  '  -  '  '  I  -  vr’  '  00 

^2v(w*  -^Xw  ♦  « cot  8)  +  PH(^  - +  cot>  8(u  -  +  2  vcot »(«  -  -  ^)!1  sin  6 


00  00'  00 
+2«i’p,H'sin  8, 


00  00  00' 


with  corresponding  differential  equations 


0  = 


0L,  0  8L|  0  0L, 

0U  00  0Uf  dt  0U, 


and 


8L,  0  8L,  0  8L, 

0W  00  0tVf  0t  0W, 


(17) 

(IB) 

(19) 


Substituting,  we  find 
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0  =  (U^’)  + 

+((1  +  V)  +  ^’(v+  co»*  ~ 

Eh  36'  36'  Se 

-  0^cotd(2-  v'  +  cot*fl)-^-2(l+  v)w-—(l  +  fi')w. 


Differential  equations  (20)  and  (21)  have  solutions  of  the  fonn 


u(»?)=Xl/.(l-»l‘)'"^. 

••0 


«'(»!)=  in' /’.(n).  (23) 

••0 

where  ri  =  cos6  and  E,(r})  are  the  Legendre  polynomials  of  the  first  kind  of  order  n.  The  differential  equations  of 
motion  (20)  and  (21)  are  satisfied  if  the  expansion  coefficients  U,  and  M'  satisfy  a  homogeneous  system  of  linear 
equations. 

Vacuum  Case 

If  we  consider  the  simpler  vacuum  case  first,  where  =  0,  the  linear  equations  are 

0  =  in’(l  +  6^')-(l  +  ^*)icK/.  -[^'()r-3n’)+l+  v]W„  (24) 

ad 

O  =  -2.09’ir  +  l+  vX/, +(n’(l  +  ^*)-2(l  +  V')-^’»rA.)H'..  (25) 

where  n  =  wo /c,.  r  =  v  +  A.  - 1.  and  A.  =  «(n  + 1).  The  determinant  of  (24)  and  (25)  yields  a  frequency  equauon 


of  the  form 


0  -  Q^(l  +  7^')  +  n'Hl  +  2^*)r  - 12(1  +  vXf  +  6^*)  +  P'kX,] 
+3^'A,  (1  ■»•  V))  +  (A.  -  2X1  -  1^)  ♦  PM«rA.  ■*•  2(1  -  A.Xl  +  v)l 


We  begin  consideration  of  the  fluid  loaded  case  by  noting  that  "for  a  plate,  fluid  loaded  on  one  side,  of  mass  per  unit 
area  p,h,  the  appropriate  nondimensionai  measure  of  fluid  loading  at  a  frequency  to  is  pc  /  <op,h.  (lunger  and  Feit, 
1986,  p.  237).  Analogously,  we  may  expand  the  surface  pressure  for  a  sphere  in  terms  of  modal  specific  acoustic 
impedances  r,  as  follows. 


p(a.  fl.  #)=  X  («K  ®)ca  m#. 


.  h.{ta) 

z,  -tpc— - . 

Kita) 


Splitting  z.  into  real  and  imaginary  parts,  we  have 


».  =  r.  -  •‘i*".. 


1 '>;(*-)] 

and 

w  I  *,'(*«)  J 

For  our  simpler  case  of  nontorsional  cnsonificatioo.  the  surface  pressure  expansion  simplifies  to 
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^>.W  =  -'^t,w,p,{ccse), 

m»0 


which  by  substitution  becomes 

PtW  =  -  ta’tVm, )P,(cos0). 

<i>0 

Substitution  of  (33)  into  (20)  and  (21)  will  result  in  simultaneous  linear  equations  of  the  form; 

0  =  (QMl  +  6^’ )  -  (1  +  )r -  US*  (If  -  3ft' )  f  1  +  v)W. . 


and 


(32) 

(33) 

(34) 


0  =  -A,,(^*if  +  l+  v)t/. +(ft*(l  +  -^^^  +  ^*)  +  i--^n-2(l+  v)-^**-A.li3'..  (35) 

p,h  k  p,c. 


Setting  the  real  part  of  the  detenninant  of  (34)  and  (35)  to  zero  results  in  a  quadratic  equation  in  ft'for  the  fluid 
loaded  case.  Ifwedeflne 


a  =■ 


(36) 


artd  neglect  terms  of  order  greater  than  then  the  quadratic  is,  rmally, 

0-ft*(l  +  o  +  7^*  +6a^*)-ft*((l  +  a  +  2^*  +  a/3*)«-  +  2(1+  v)(I  +  6)?*)  +  -  3^*2.(1  +  v)) 

+3^*A.(1  +  v))  +  (A.  -  2X1  -  s^)  +  ^’«r(rA.  +  2(1  -  A.Xl  +  v)J.  (37) 


CX)NCLUSIONS 

With  (26)  and  (37)  in  hand,  the  obvious  next  step  is  to  test  the  results  numerically  against  exact  results.  We  note 
that  by  setting  to  zero  we  revert  our  solutions  to  previously  derived  models  (lunger  and  Feit,  1986),  S imihrly, 
setting  a  to  zero  in  (37)  is  equivalent  to  removing  the  fluid  loading  from  the  model.  This  makes  (37)  revert  w 
(26).  By  altentately  retaining  or  zeroing  /S’  in  (26),  and  similarly  for  a  in  (37),  we  have  three  distinct  models  with 
differing  degrees  of  physicality.  We  can  use  the  resulting  expressions  to  calculate  resonances  over  a  frequency  range 
and  compart  them  with  the  exact  results.  We  may  then  rank  the  various  af^ximaiions  in  order  of  their  agreement 
to  the  exact  results.  Limiiauons  of  each  of  the  methods  can  then  be  outlined  as  well  as  those  of  sheii  methods  m 
general. 
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